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ABSTRACT 


The accuracy of user navigation fix based on the NAVSTAR 
Global Positioning System is described by a 4x4 position- 
time error covariance matrix. The "trace" of this matrix 
serves as a convenient navigation performance index and 
the square-root of the trace is called Geometric Dilution 
of Precision (GDOP) . In this paper, certain theoretical 
results concerning the general properties of the navigation 
performance are derived. An efficient algorithm for the 
computation of GDOP is given. Applications of the results 
are illustrated by numerical examples. 
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INTRODUCTION 

The NAVSTAR Global Positioning System (GPS) , when fully op- 
erational in the early 1990's, will provide world-wide 
navigation through synchronized transmissions from a con- 
stellation of eighteen 12-hour period satellites in three 
55°-inclination orbital planes. An accurate user navigation 
fix (position and time) will be obtainable by receiving 
transmissions from four satellites and decoding the signal 
transit times. 

One may relate the measurements, referred to as the pseudo- 
ranges, to the navigation state as follows 

°Tj = Jcx ,-*.) 1 + (X*~*,y* + X* + (i) 

where C - velocity of light 

Tj = Signal transit time from GPS satellite 
"j" to user, not corrected for user 
clock offset, At 

X. X 0 X_ X. = user naviagation state, the first three 
represent a set of convenient Carr.esian 
user coordinates, X^ = CAt is a range 
bias equivalent of user clock offset 
x 2 x^ = corresponding Cartesian coordinates of 
GPS satellite "j" 
n^ = random measurement noise 

— 

Senior Principal Engineer, Orbit Operations, System 
Sciences Division. 
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From a set of four measurements, a user navigation fix may 
be determined. The accuracy of the fix is characterized by 
the following 4x4 position-time navigation error covariance 
matrix 
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(H WH) 
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a, b, c and d = line-of-sight unit vectors from a set 
of four GPS satellites to the user, W = 4x4 covariance 
matrix of random measurement noise, superscript "T" = 
transpose of matrix. 

The measurement error covariance matrix W is generally taken 
to be diagonal, which is strictly true for uncorrelated 
measurements only. In practice, assignments of quantitative 
values to the elements of W also takes into consideration 
such factors as the elevation and health status of individual 
GPS satellites. Thus W may be more appropriately be re- 
ferred to as the weighting matrix. For uniform weighting, P 

T — 1 

is proportional to (H H) , which depends only on the rela- 
tive geometry of the user and the four GPS satellites, as is 
evident from Equation ( 3 ) • The square-root of the "trace" 
of (H H) is referred to as Geometric Dilution of Precision 
(GDOP) , a self-explanatory name. Whatever the weighting 
strategy, the "trace" of the navigation error covariance 
matrix serves as a covenient and natural performance index 
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characterizing the accuracy of the naviagation fix. For 
a diagonal weighting matrix W, 

T -1 

TRACE "P" = sum of diagonal terms of (H H) weighted 
by the inverse* of the corresponding ele- 
ments of W 

Thus the evaluation of the GPS naviagation performance is 

essentially equivalent to the computation of the diagonal 
T -1 

terms of (H H) , which may be called the GDOP matrix for 
convenience . 

The navigation performance index. Trace "P" , also serves as 
a criterion for the selection of a set of four best GPS 
satellites among those visible, which may be as many as ten 
for users which are satellites themselves. If, for optimum 
performance, each of the different combinations of four has 
to be evaluated, the computational burden can be considerable. 
In the following, certain theoretical results concerning 
the general properties of the GDOP matrix are derived. An 
efficient algorithm for the computation of GDOP matrix and 
the navigation performance index is given. Applications of 
the results are illustrated by numerical examples. 

ANALYTICAL RESULTS 


To solve for a navigation fix from four measurements, the 
partial derivative matrix H must be non-singular. Since 


determinant H = 


a T -d T 

b T -d T 

c T -d T 


0 

0 

0 

1 


a -d 
b -d 
c -d 


12-4 



a navigation fix can be determined from four GPS satellites 
with line-of-sight directions a, b, c, d, if and only if 
the three vectors (a-d) , (B-d) , and (c-d) are linearly inde- 

pendent, i.e., non-coplanar. This shall be assumed to be 
the case in the following development. 

Since Trace (H^H) ^ = Trace (HH^) ^ , by making use of 
Equation (3) and the fact that a, b, c, d, are unit vectors, 
one obtains. 


Trace (H T H) 1 = Trace (HH T ) 1 


= r, 


race 
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The advantages 
become obvious 


T T 

of dealing with HH instead of H H will 
below. 


The following may be observed from Equation (4) : 

T 

1. The matrix HH in Equation (4) is non-negative, 

symmetric, and with identical diagonal terms which 

are greater than the off-diagonal terms. (Expres- 
T 

sions such as a b are scalar product of unit 
vectors and are less than unity) . These properties 
give rise to good behaviour in numerical operations. 
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Since the Trace of a matrix is equal to the sum 
of its eigenvalues and the eigenvalues of the 
matrix inverse are inverses of the eigenvalues of 
the matrix itself, one has the following results: 

a. Trace (HH T ) -1 = - + ^ ^ + jf' where the 

\L *2 *3 *4 

T 

Vs are eigenvalues of (HH ) with 

\l + V‘ A 3 + \ = 8 - 

b. From "a" above and the fact that the )^s are 
non-negative, one may conclude that 

Tmce (nk t )-'j2 < 5 

T 

c. Let us order the eigenvalues of HH as 

o < $ A x ^ X 3 

One has the obvious inequality 

True* ( M H ) 1 - \ ( ^ ( 


K f 


/ OtoCr btund + f ) 

v Xj Aj A (f. > 

Aj + Au. ^ $ 



Thus knowledge of the smallest eigenvalue of 
T 

HH provides another lower bound for the 
navigation performence index. Sometimes this 
lower bound also serves as a good estimate. 


d. The 2x2 principal submatrix of HH , e.g. 
2 


T 

a b+1 


|a T b+l 
l-a T b. 


has eigenvalues 3+a b and 


From the Theorem of Root Separation 

for Symmetric Matrices'*' one obtains the follow- 

T 

ing bounds on the eigenvalues of HH 


^ £ 1-a b s X ^ 


(7) 


\ 3 3 + a b£ V. 


( 8 ) 


These inequalities have no preferences for the 
labeling of the unit vectors. That is, a, b 
may be replaced by c, d, etc., to obtain sharper 
bounds. In particular, one must have X, < J 
and ^ • Therefore, the eigenvalues of 

HH^ cannot be all identical and the equality 
sign in (5) may be deleted. Physically, this 
follow from the fact that the four unit vectors 
in three-dimensional space cannot play identical 
roles in the four-dimensional position-time 
space. Combining inequalities (6) and (7) , one 
obtains another inequality. 


Trace (HH T ) ] > J + ^coi© (9) 

where fr = smallest angle subtended by two 
line-of-sight vectors. 


12-7 


This inequality, although not sharper than 
Inequality (8), is easier to calculate, and 
expresses the intuitive rule of thumb that an 
accurate navigation should not rely on a GPS 
constellation that is clustered together. We 
shall see later that with good geometry, navi- 
gation performance index of magnitude less than 
3 may be obtained. On the other hand, as indi- 
cated by Inequality (9) , a navigation perform- 
ance index in excess of 8.5 would result if any 
two line-of-sight vectors to GPS satellites are 
separated by 30° or less. 

An upper bound for the navigation performance 
index may be obtained 


W* C Kff T ) H x, + i {10) 

6y A, $ \ x $ Aj A* > . 

It may also be pointed out that because the determinant of 
a matrix product is equal to the product of the individual 
determinants, and that the determinant of a matrix is equal 
to the product of its eigenvalues, one has the relation 


— ( dzftr/i\r hdtd M )* — 


(t •*» <f 
fc «- cf 
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The maximization of 


has been suggested as a 


,C~cf. 3 

convenient GPS selection criterion . It is seen from the 
above equation that this criterion is equivalent to a maxi- 
mization of the denominator of our performance index. 
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ALGORITHM 


An efficient algorithm for the computation of the GDOP matrix 
may be obtained by noting the following decomposition of the 
measurement partial derivative matrix: 
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From this decomposition, the Sherman-Morrison Formula gives 
us 

A' 1 - 


(ID 


Then one has, by straight-forward simple algebra, 


A = 


f ! t\ ^ ! -*T 


° I ' 


and 


H = 


oi 


~d T f ' -d\ ! 


* 


(13) 


where q = f + g +h 
cC — 1 - d T q, 

When H -1 is obtained, one may obtain the GDOP matrix as 
(HH T ) -1 = (H -1 ) T (H -1 ). In particular, 

m — 1 **4 

Trace (HH ) = sum of the squares of the elements of H 

= f ' f + ff +h T h + i { C/f) ft 4 + (^)fhj 

+ { -+ (d r f )‘+ (<•>)*+ f j t a)*J < 14 > 


Equations (11), (13) and (14) constitute the algorithm. It 

reduces the reversion of the 4x4 matrix HH to the inversion 
of a 3x3 matrix (a [ b ! c) plus the scalar products of sev- 
eral 3x1 vectors. Notice that Eq. (13) may also be obtained 

2 

from inverting K by partitioning . But the Sherman-Mornson 
Formula provides additional flexibility as will be discussed 
below . 
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An important advantage of this algorithm is that very little 
recomputation is required when the fourth GPS satellite is 
switched. In selecting the best set of four GPS satellites 
from the many possible combinations, a simple combinatorial 
test logic may be advantageous. For this purpose, one may 
need the flexibility of changing any one of the rows of H. 
Although Eq. (14) remains valid provided one interprets the 
vectors f, g, and h accordingly, this does mean these vec- 
tors have to be recomputed. In that case it is preferable 
to use Eq. (11) directly instead of Eq. (14) . To illustrate 
let us assume that for a particular GPS configuration, 

H _ ^=G is already obtained. If the nth (n = 1, 2, 3, 4) GPS 
Satellite with line of sight vector r is to be replaced by 
another satellite with line-of-sight vector p, the new 
measurement partial derivative matrix may be written as 


H = H + 


<5V 


30 


& 


4n 


(p-r) T j o 


where <5^ is the Kronecker delta ( f° r n $in ^ ^ or 


i = n) . 


From the above decomposition the Sherman-Morrison Formula 
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The computational economy provided by this equation is ob 
vious . 
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APPLICATIONS 


ihbitively, orthogonal line-of-sight user-to-GPS satellite 
configurations are favorable. In three-dimensional space, 
it is, of course, impossible to have a set of four mutually 
orthogonal unit vectors. An alternative has three of the 
line-of-sight vectors a, b, c orthogonal. For this case 
the vectors f, g and h become the same orthogonal unit vec- 
tors as a, b, c and Eg. (14) simplifies to 


T r ac* rNK T r'= 3 t A ( d, + ^ + d 3 )/<* t * ( M- + J/ )/ u * 

— 1 + { $ + 2 

2 2 2 

where d^,d 2 fd^ with d^+d 2 +d^ = 1 are components of the line- 
of-sight unit vector d along the orthogonal a,.b, c direc- 
tions. ...It is of interest to note that for this case the 
Navigation Peformance Index depends only on (d.+d-td.), the 

_L ^ J 

simplest symmetric function of the components of the vec- 
tor d. The best performance index of 2.80 is achieved for 
T 

d = (-lf~l,-l)//T. This is the situation that the line of 

sight to GPS satellite "d ,! shows no preference to, but is 

directed away from the other GPS satellites, an artificial 

but not improbable configuration for an user satellite. 

T 

For d = (lfl,l)/x/3, i.e., d having the same general direction 
as the other three lines-of-sight , the performance index 
degrades to 13.20. This degradation reminds us of the state- 
ment made earlier about avoiding closely-grouped GPS satel- 
lites. For d = -a, i.e. for an user located between two 
GPS satellites, the performance index has the value 4.00. 

There is reason to think that a GPS constellation with a-d, 
b-d, c-d orthogonal may give good navigation performance. 
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This may be realized with the set of line-of-sight vectors 

*7 * i/O * T = ( \, -i , • ) /Js t c T ~ Co »,-») A/J 

and d T = (1, 1, l)/,/7. However, for this configuration, 
the angle between the vector d and any other vector is 
cos - '*' (2/3) , which is comparatively small, and may be undes- 
irable from the consideration of the preceeding section. 
In4e<?J , it follows immediately from Inequality (10) that 
the navigation performance index must be in excess of 
9/8 + 1/(1 - 2/3) = 4 '/% , a lower bound which may be com- 
pared with the exact index of 5.5 obtainable from straight- 
forward simple computation. On the other hand, by reversing 
the direction of the vector d given above, one has the com- 
pletely symmetrical configuration that the line-of-sight 
vectors are all separated by the same angle cos (- 1/3) . 

For this configuration one may compute the eigenvalues of 
HH T as \ 1 = K 2 = \ 3=4/3 and \ 4 = 4, giving rise 
to a navigation performance index 

y -t . A . _L _ a x — 

(race (MU') ’ *= * * J T > r V 

Notice that for this configuration, 

1. The upper bound for X, given in Inequality (7) 
is achieved. 

2. Any perturbation of the configuration will result 
in a decrease in the minimum angle between two 
line-of-sight vectors , and therefore a decrease 
in \ x . 
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Thus this configuration maximizes the smallest eigenvalue 
T 

of HH , or equivalently, minimizes the largest eigenvalue 
T -1 

of (HH ) . Whether this also happens to be the best con- 

figuration remains to be investigated. 
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